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Abstract 

This paper is concerned with open quantum systems whose dynamic variables satisfy canonical commutation 
relations and are governed by quantum stochastic differential equations. The latter are driven by quantum Wiener 
processes which represent external boson fields. The system-field coupling operators are linear functions of 
the system variables. The Hamiltonian consists of a nominal quadratic function of the system variables and 
an uncertain perturbation which is represented in a Weyl quantization form. Assuming that the nominal linear 
quantum system is stable, we develop sufficient conditions on the perturbation of the Hamiltonian which guarantee 
robust mean square stability of the perturbed system. Examples are given to illustrate these results for a class 
of Hamiltonian perturbations in the form of trigonometric polynomials of the system variables. 

Index Terms 
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I. INTRODUCTION 

The quantum mechanical concept of quantization is concerned with assigning quantum observables to 
classical variables (and functions thereof). Weyl’s proposal for the development of a general quantization 
scheme was introduced in 1927 (see for example, [23, Section IV. 14]) soon after the invention of quantum 
mechanics. An important feature of the Weyl association is that it treats quantum dynamic variables 
equally and leads to correct marginal distributions for them [1, Chapter 8]. The Weyl quantization 
scheme employs Fourier transforms and is known to be a convenient and, in many respects, satisfactory 
procedure for quantization [1], [3], [24]. In addition to providing a mathematical formalism, this scheme 
also offers an interpretation of quantum mechanical phenomena, thus leading to a better understanding 
of their physical aspects [1], [24]. The aim of the present paper is to use the Weyl quantization for the 
modeling of perturbations of Hamiltonians for a class of open quantum systems and the robust stability 
analysis based on this description of uncertainty. 

A wide range of open quantum systems, which interact with their environment, can be modelled by 
using the apparatus of quantum stochastic differential equations (QSDEs) [7], [14]. In this framework, 
which follows the Heisenberg picture of quantum dynamics, a quantum noise is introduced in order 
to represent the surroundings as a heat bath of external fields acting on a boson Fock space [14]. The 
QSDE approach to open quantum systems is employed by the quantum dissipative systems theory [8] 
which addresses robust stability issues. 

The robustness of various classes of perturbed open quantum systems, modelled by QSDEs, has been 
addressed in the literature using dissipativity theory and different notions of stability (see for example 
[15], [16], [20], [21]). In particular, robust mean square stability with respect to a class of perturbations 
of Hamiltonians has been studied in [15] and its applications have been presented in [17], [18]. In these 
papers, the classical and quantum models of the perturbed Hamiltonian are related by using power series 
of quantum variables in combination with Wick’s quantization [1, pp. 445]. 

In the present paper, we consider a class of open quantum systems whose dynamic variables satisfy 
Heisenberg canonical commutation relations (CCRs) and are governed by QSDEs. The system-field 
coupling operators are assumed to be known linear functions of the system variables, while the 
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Hamiltonian is split into a nominal quadratic part and an uncertain perturbation which is, in general, 
a non-quadratic function of the system variables. Following a similar approach in [22], we use Weyl 
quantization in order to model the perturbation of the system Hamiltonian. The fact that the Weyl 
quantization employs Fourier transforms makes it particularly suitable for modelling uncertainties in the 
form of trigonometric polynomials of the system variables, such as in [18]. Assuming that the nominal 
linear quantum stochastic system is stable, we develop sufficient conditions for the robust mean square 
stability of the perturbed system. These conditions employ a linear matrix inequality (LMI) which, 
in addition to the conventional Lyapunov part, involves a linear operator acting on matrices, whose 
structure is analogous to that of generalized Sylvester equations (see, for example, [4]). Such operators 
play an important role for moment stability of quasilinear quantum stochastic systems [21]. 

The rest of the paper is organized as follows. Section Hill describes the class of open quantum systems 
being considered. Section [IV] models Hamiltonian perturbations in a Weyl quantization form. Sections lYl 
and |Vj discuss the time evolution of weighted mean square functionals of the system variables and a 
related dissipation inequality. Section IVIII provides sufficient conditions of robust mean square stability 
for an admissible set of Hamiltonian perturbations. Section IVIIII discusses techniques for verifying 
these conditions in terms of the Weyl quantization model. Section HXl provides examples to demonstrate 
applicability of the approach. Section |X] summarizes the results of the paper. 

II. NOTATION 

Unless specified otherwise, vectors are organized as columns, and the transpose (-) T acts on matrices 
with operator-valued entries as if the latter were scalars. For a vector X of operators X\,...,X r and a 
vector Y of operators T|,..., Y s , the commutator matrix is defined as an (r x v)-matrix [X,Y T ] :=XY T - 
(TA t ) t whose (j,k )th entry is the commutator [Xj, p] XjY^ — YkXj of the operators Xj and p. Also, 
(•)"' := ((-) # ) T denotes the transpose of the entry-wise operator adjoint (-) # . In application to complex 
matrices, (-) 1 reduces to the complex conjugate transpose (•)* := ((-)) T . Furthermore, S r , A r and H r := 

+ iA r denote the subspaces of real symmetric, real antisymmetric and complex Hermitian matrices of 
order r, respectively, with i \/—T the imaginary unit. Also, /, denotes the identity matrix of order r, 
positive (semi-) definiteness of matrices is denoted by (^) >-, and <g) is the tensor product of spaces or 
operators (for example, the Kronecker product of matrices). The adjoints and self-adjointness of linear 
operators acting on matrices is understood in the sense of the Frobenius inner product ( M,N) Tr (M*N) 
of real or complex matrices, with the corresponding Frobenius norm ||M|| := \J(M.M) which reduces 
to the standard Euclidean norm | • | for vectors. Also, 11 v | p := Vv T Kv denotes the Euclidean norm of a 
real vector v associated with a real positive definite symmetric matrix K. Finally, E<§ := Yr(pq) denotes 
the quantum expectation of a quantum variable q (or a matrix of such variables) over a density operator 
p which specifies the underlying quantum state. For matrices of quantum variables, the expectation is 
evaluated entry-wise. 


III. OPEN QUANTUM STOCHASTIC SYSTEMS 

We consider an open quantum stochastic system interacting with an external boson field. The system 
has n dynamic variables X\,...,X n which satisfy CCRs 


[X,X T ] = 2/0, 


PCI 


(l) 


where the CCR matrix 0 e A„ is assumed to be non-singular. The system variables evolve in time 
according to a QSDE 


dX 


i[H,X] - 1 A j dt + BdW. 


( 2 ) 
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Here, W is an m-dimensional vector of quantum Wiener processes W\ ,..., W, n with a positive semi- 
definite Ito matrix Q. G H m : 

dWW T = Qdt, Q. = I m + iJ , (3) 

where J e A m . The matrix B e R" xm in ([2]) is related to a matrix M e R" 7X " of linear dependence of 
the system-field coupling operators on the system variables by B — 2QM J . The term 
in the drift of the QSDE £1]) is the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) decoherence 
superoperator [6], [11] which acts on the system variables and is associated with the system-field 
interaction. Also, H is the Hamiltonian which describes the self-energy of the system and is usually 
represented as a function of the system variables. For what follows, we assume that H is split into two 
parts: 

H = H 0 +Hi. (4) 


Here, 

H 0 := ] -X J RX = i £ r jk XjX k (5) 

2 Z j,k= 1 

is a quadratic function of the system variables with a real symmetric matrix R {rjk)\< tJ ,k<n °f order 
77 , which corresponds to a nominal open quantum harmonic oscillator [2], [5]. Also, H\ is a self-adjoint 
operator on the underlying Hilbert space, which is interpreted as a perturbation of the Hamiltonian and 
is described in Section [TV] By substituting © and ([5]) into © and using the CCRs ©, it follows that 
the QSDE takes the form 


dX = {AX + Z)dt + BdW, Z := i[H h X] , (6) 

where Z is an /7-dimensional vector of self-adjoint operators, and A e R' !X " is given by 

A:=2@R--BJB t Q- 1 . ( 7 ) 

2 

It is assumed that the matrix A is Hurwitz, and hence, the nominal open quantum harmonic oscillator is 
stable. In particular, in the absence of perturbations (that is, when H i = 0, and the system is governed 
by a linear QSDE dX —AXdt+BdW), the system variables have finite steady-state moments of arbitrary 
order. Also, note that the CCRs (Q]) are preserved in time for any perturbation of the system Hamiltonian. 
This is a consequence of the joint unitary evolution of the dynamic variables of the system and its 
environment, with the CCR preservation being part of the quantum physical realizability conditions [9]. 

IV. PERTURBATION OF HAMILTONIAN 

We will model the perturbation Hi of the system Hamiltonian in © by using the Weyl quantization 
[3] as follows: 

Hi := [ h{Z)e xTx dZ. ( 8 ) 

J R» 

Here, h : M' ! —* C is a complex-valued function which satisfies h{—X) — h{X) for all A £ R'\ thus 
ensuring that H\ is a self-adjoint operator on the underlying Hilbert space. The function h is computed 
as the standard Fourier transform 

h{X) := (27r)“ n [ Hi{x)e~ aTx dx (9) 

J R” 

of a real-valued function H\ : R' 1 —» R of n classical variables whose quantization leads to ®. The 
following lemma computes the perturbation term in the governing QSDE in the presence of perturbations. 
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Lemma 1: The perturbation vector Z in the drift of the QSDE ©, which corresponds to ([8]), can be 
computed as 

Z = 2/0 f h(X)Xe xTx dA. (10) 

JM. n 

Proof: By substituting ([8]) into the definition of Z in © and using the bilinearity of the commutator, 


it follows that 

Z = i f h{X)[e aTx ,X]dX. 

J K” 

(ID 

Now, 

[e iiT)I ,X] = e iiTx X-Xe‘ lTx 



= (t ,llx Xt- ,llx -X)t ilTx = 2©Ae' ATx , 

(12) 


where use is made of the relation 

e aT]C Xe~ aTx =e [aTx ’\x) = £ -] fc (X) 

k=o ^ 

= X + [iX T X,X]=X-i[X,X T ]X = X + 2@X 

which follows from Hadamard’s lemma [12] and the CCRs (Q]). Here, {£,, -] k denotes the A'-fold application 
of the commutator with a given operator Substitution of (fT2l) into (fTTI) leads to (flOl) . ■ 

Note that the vector Z in (flOl) can be regarded as the Weyl quantization of the M' 7 -valued function 
jc i— 2®d K H\ (x), where d x (-) denotes the gradient operator. Indeed, from © it follows that ih(X)X is 
the Fourier transform of the function d x H\. 

V. WEIGHTED MEAN SQUARE FUNCTIONALS 
For what follows, consider a weighted mean square of the system variables: 

V := E(A T IIA) = £ 7C jk E(XjX k ), (13) 

j,k= 1 

where n := (Jtjk)i^j,k^n £ is a given matrix. Here, the quantum expectation Ec, := Tr(p^) of a 
quantum variable t, is taken over a density operator p with a product structure p — U5®v, where G7 is 
the initial quantum state of the system, and v is the vacuum state of the external boson fields. We will 
now consider the time evolution of the quantity V. 

Lemma 2: For the perturbed quantum system governed by ©, the quantity V in (fT3T) satisfies a 
differential equation 

V =<A T n + nA,/ 5 ) + <n 1 RR T ) 

+ <n,ReE(AZ T + ZA T )>, (14) 

where 

P:=ReE(XX T ) (15) 

is the real part of the matrix of second moments of the system variables. 

Proof: Consider a matrix S 6 EI„ of second moments of the system variables given by 

S:=E(XX t ) =P + z'0, (16) 

where use is made of (|T|) and (fl5l) . Since n in (fl3l) is a real symmetric matrix, then 

v = (n,s) = (n,p). (17) 
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By combining the quantum Ito lemma and the Ito product rules with ©, it follows that 

d(ZX T ) = (dA)X T + AdX T + (dA)dA T 
=((AX+Z)dt + £dfF)X T 
+ X ( (X t A t + Z T )df + dl¥ T B T ) 

+ BdWdW T B T 

=(AXX t +XX t A t +BQB t + XZ t + ZX T )dt 
+ BdWX T + XdW T B T , (18) 

where use is also made of ©. Since the external boson fields are assumed to be in the vacuum state, 
the averaging of the QSDE (fT8T) leads to a differential equation for the matrix S in (fT6l) : 

S = AS + SA t +BQ.B t + E(XZ t +ZX t ). (19) 

In view of ©, (fl6l) and (fT71) . the substitution of (fl9l) into V = (II,ReS) and using the duality relation 
(n,AP + PA T ) = (A T n + nA,P) leads to (Q3J, thus completing the proof of the lemma. ■ 

Note that the first line in (fl4l) corresponds to linear dynamics of the nominal open quantum harmonic 
oscillator, while the second line comes from the perturbation of the system Hamiltonian. 

VI. DISSIPATION INEQUALITY 

Following [20, Section VIII], we say that a matrix L (Lj k ) | < h k<:r of linear operators on the 
underlying Hilbert space, satisfying L T = L, is superpositive if the self-adjoint operator u*Lu \— 
YJj k= \ TTjU k Ljk is positive semi-definite for any vector a («j)t^ r £ C r . This notion extends the 
standard positive semi-definiteness from a single operator to a matrix of operators and will be written 
using the same symbol as L ip 0. Note that the superpositiveness Lip 0 implies that EL ip 0 for any 
density operator over which this expectation is taken. Indeed, Lip 0 implies that uELu = E (u*Lu) ^ 0 for 
any u £ C r . Similarly to the usual positive semi-definiteness for single operators, the superpositiveness 
induces a partial ordering for matrices of operators. An example of a superpositive matrix is provided 
by XX T because u*XX T u — ( u*X)(u*Xy ip 0 for any u £ C". Now, consider an operator inequality 

d 

zz T ^ Ml £ nxx T r T k + juo I n (20) 

k= 1 

for the vector Z from © in the sense of superpositiveness. Here, r|..... £ M' !X " are fixed matrices, 

Pi and po are real constants, and d is a positive integer. The following lemma derives a useful upper 
bound from (l20l) . 

Lemma 3: Suppose the vector Z satisfies (l20l) . with p\ > 0. Then 

xz T +zx T ^ £ r k xx T rJ + — /„, r 0 := ^p[i n . 

k =0 

Proof: Consider an auxiliary vector Y of self-adjoint operators defined by Y ^fpfX 
Since the matrix of operators YY T = p\XX T +■ jpXZ T —XZ T —ZX T is superpositive, then 

T T T 1 T 

XZ T + ZX T =4 piXX T -i -ZZ T 

Ml 

^ PiXX T + £ r k xx T r] + ^i n . 

k= i 

The last inequality follows from (l20l) and leads to (12TI) . 


( 21 ) 
1 z 

v/Fi 
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We will now use Lemma [3] in order to obtain a dissipation inequality for the mean square functional 
of the perturbed system. 

Lemma 4: Suppose the vector Z in © satisfies the operator inequality (l20l) . with /ii > 0. Also, let 
the weighting matrix n in (fl3l) be positive semi-definite and satisfy the following LMI 

d 

A T n+m + £rJnr fc +yn^o, (22) 

k =0 

where y is a real constant, and To is the matrix given by (I2TT) . Then the quantity V in (fl3l) satisfies a 
dissipation inequality 

V < -yV + (n,BB T ) + — Trn. (23) 

Mi 

Proof: From ([22]) and the positive semi-definiteness of the matrix P in (fT6l) . it follows that the 
first inner product in fli~4l) admits an upper bound 


<A T n+n4,p) < 


rjnr, 


yn,P 


k =0 


(24) 


By applying Lemma [3] and using the monotonicity of the quantum expectation with respect to the 
superpositiveness, it follows from (l2TT) that 


E(XZ T +za t ) ^ £ E(r k xx T rl) + 

k=o Mi 



Here, use is made of (fl6l) and the fact that To,.... Y ( i are constant matrices. Therefore, since H 0, 
then (1251) leads to the following upper bound for the last inner product in 


<n,R e E(« T + z* T )><(n,£r^ + ^) 

LrJnr,,p\ + ^Trn. (26) 

\k=0 / Ml 

It now remains to note that substitution of (l24l) and (l26l) into (Ti~4l) establishes (1231) . ■ 

Note that, in addition to the usual Lyapunov part A T n-t-IL4 (which comes from the nominal linear 
system), the LMI (l22l) involves a linear operator 

d 

n^^rjnr, (27) 

k =0 

on the space S„, which, in view of (l20l) . is associated with the perturbation of the Hamiltonian. The 
structure of the linear operator in (l27l) is analogous to that of the generalized Sylvester equations [4]. 
Such operators are present in moment stability conditions for quasilinear quantum stochastic systems 
[21, Section IX]. Furthermore, this operator structure resembles the Kraus form of quantum operations 
[13, pp. 360-373], 
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VII. MEAN SQUARE STABILITY 

The following lemma provides sufficient conditions of mean square stability [15], [16] of the perturbed 
quantum system ©. 

Lemma 5: Suppose the vector Z in ([ 6 ]) satisfies the operator inequality (l20l) for some /ii > 0, 
ri,...,r d G W ixn and Mo G R. Also, suppose there exist a weighting matrix 11 >- 0 and a constant 
y > 0 which satisfy the LMI (1221) . Then the quantum stochastic system ([ 6 ]) is mean square stable, with 
the upper limit of the quantity V in (fT3l) satisfying 

limsupV ^ -( (T1,BB t ) + —Trll). (28) 

yv x Mi > 

Proof: The upper bound (l28l) follows from the Gronwall-Bellman inequality applied to (1231) : 

V(t) f V(0)e“^ + ( <n,BB T ) + ^TrlTj ~~~~~> 

which holds for all times t ^ 0, where use is made of the integral f 0 c~ r 'd.v = -—y—. It now remains 
to note that ( 1281 ) implies mean square stability of the quantum system being considered in view of the 
assumption that II >~ 0, whereby E(A T A) < i~ 7 ny’ with A m j n (fl) > 0 denoting the smallest eigenvalue 
of the weighting matrix. ■ 

In combination with the assumption of stability for the nominal linear quantum system, Lemma [5] 
leads to the following sufficient conditions on the perturbation Hamiltonians which guarantee mean 
square stability of the perturbed system. 

Theorem 1: Suppose the matrix A in ([7]) is Hurwitz. Then the following set is nonempty: 

^:={([ii,r 1 ,..,r rf )Gix(r x Y: d>o, m >o, 

(122|) is satisfied for some n >~ 0, y > oj. (29) 

Furthermore, the perturbed quantum system ([ 6 ]) is mean square stable for any perturbation Hamiltonian 
H i from the uncertainty set 

:= j//| given by ([ 8 ]) : Z satisfies (1201) 

for some (pi,r u ...,r d ) G &, /io G m|, (30) 

where Z is the vector of operators associated with H\ by ©. 

Proof: Since the matrix A is Hurwitz, then there exist II 0 and 70 > 0 such that A T n + TL4 + yon ^ 
0. Hence, by continuity, there exist (for any given positive integer d) sufficiently small > 0 and 
r h ...,T d G M. nxn such that the LMI (l22l) is satisfied for the same n and a smaller y > 0 (that is, 
0 < y < yo). This proves that the set f J ? in (l29l) is indeed nonempty. Now, the property, that the perturbed 
system (| 6 ]) is mean square stable for any perturbation Hamiltonian H\ belonging to the class in (l30l) . 
was established in Lemma [5} ■ 

VIII. TECHNIQUES FOR VERIFYING THE OPERATOR INEQUALITY 

We will now discuss several techniques for verifying the operator inequality ([20]) which plays a central 
role for the uncertainty class in (l30l) . The following lemma reformulates this operator inequality in 
terms of the Weyl quantization of the perturbation Hamiltonian. 

Lemma 6: The operator inequality (l20l) for the vector Z in © is representable in the form 

4©/ fi(/l)e aT *dA0 ^ Ml Yr k XX T r T k +HoI n . (31) 

J R" 
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( 32 ) 


Here, the function h : M" —> C nxn is expressed in terms of the Fourier transform h from ([9]) as 

//(A) := [ //(A - t)//(t)(A - z)T T e nT@?L dT. 

_ . J w 

Proof: From (1101) and the antisymmetry of the CCR matrix 0, it follows that 

ZZ T = 40 f h{ A)Ae /AT *dA / //(rjrV^dT© 

,/R" J R" 

= 40 / h{k)h{T)XT T e aTx e itTx d?idT® 

J R 2 " 

= 40 / h(?i)h(T)?iZ T e izTm e KX+T)Tx d?idT@ 

JR 2 " 

= 40 [ //(A)e' AT ^dA 0, (33) 

JR" 

where the function h is given by (l32l) . Here, use is made of the Baker-Campbell-Hausdorff formula [12, 
pp. 40] and the CCRs (Q]), whereby 

& a T x Q iz T x = Q \\a T x,iz T x\ Q i{x+z) T x = f jz T @x ( j(x+z) T x 

and the standard change of variables (A,t) (-»• (A — T. t) for convolution integrals. Substitution of (1331) 
into (l20l) leads to (I3T1) . ■ 

The following lemma, which is given here for completeness, extends the ordering of real-valued 
functions of a real variable to the case when they are evaluated at a self-adjoint operator. 

Lemma 7: Suppose f,g : \a,b\ E arc continuous functions satisfying f(z) f g(z ) in an interval 
a ^ z ^ b. Then f(K) ^ g(K) for any self-adjoint operator K whose spectrum is contained by this 
interval. 

Proof: The spectral theorem (see, for example, [19, pp. 263]) implies that a self-adjoint operator 
K, described in the lemma, is representable as K — zv(dz), where v is a projection-valued measure on 
the interval [a,b\ C R. Therefore, since f(z) — g(z ) ^ 0 for all z G [a,b\, then the operator g(K) — f(K) = 
fa(s(z) ~ /(^)) v (d^) is positive semi-definite, and hence, f(K) ^g(K). ■ 

Note that Lemma [7] is useful for investigating the superpositiveness of a matrix of operators which 
are collinear to a real-valued function of a given self-adjoint operator. The following lemma studies a 
combined effect of several perturbations of the Hamiltonian on the operator inequality (l20l) . 

Lemma 8: Suppose the perturbation vector Z in ([6]) is decomposed as 

d 

z='£c k z k, (34) 

k= l 

where Ck are real coefficients and Z|.....Z^/ are /(-dimensional vectors of self-adjoint operators on the 
underlying Hilbert space. Also, suppose the operator inequalities 

Z&Zj =<: + IA)kIm k=l,...,d, (35) 

are satisfied for some matrices e M" x ” and real constants /ii > 0 and /iojt- Then the vector Z in (1341) 
satisfies (1201) with the same constant jl\ and the following parameters F],.... V d and /io: 

d 

:= \fo~kCk ( d > k■ Mo := ^ °k c h J -()k- (36) 

k= 1 

where the coefficients 0 / { are computed in terms of arbitrary positive scalars Vjk — Vkj as 

k— t d i 

of: = i + v jk + y; —— , a= 1 ,..., d. ( 37 ) 

j=\ j=k +1 V M 
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Proof: In view of (l34l) . 

d 

ZZ T = £ c \z k zl + £ CjC k (ZjZj +z k zj). (38) 

1 ^j<k^d 

By applying the completion-of-the-square technique, used in Lemma |3j to the vectors of operators cjZj 
and c k Z k , it follows that 

(f. 

c j c k(ZjZj + Z k zj ) ^ -^-ZjZj + VjkCiZkZj, (39) 

v jk 

where V k j — Vj k are arbitrary positive scalars. Substitution of (l39l) into (l38l) and combining the result 
with (031) leads to the following operator inequality (l20l) for the vector Z in (1341) : 

zzuiclz k zj + £ (izjZj + VjkclZkZj) 

k= 1 1 ^j<k^d V J k 

d d 

= £ o k c\z k zl ^ Ml £ r,xx T rJ+Mo4, 

k= 1 k= 1 

where use is made of the notation (l36l) and (07l) . ■ 

IX. ILLUSTRATIVE EXAMPLES 

The following examples aim to demonstrate an application of Theorem Q] to the robust mean square 
stability analysis of the open quantum system © when the perturbation Hamiltonian is a trigonometric 
polynomial of the system variables. 

Example 1: Suppose H\ cos(A ( }X), where Ao £ R” is a constant vector of spatial frequencies. This 
perturbation Hamiltonian results from the Weyl quantization ® of the function cos(A^) whose Fourier 
transform ([9]) is given by /i(A) = |(<5(A — Ao) + <5(A + Aq)), with S(-) denoting the n-dimensional Dirac 
delta-function. Substitution of h into (flOl) yields Z = —20Aosin(A ( }X), and hence, 

ZZ T = -40A o sin 2 (A o T X)A o T 0, (40) 

which can also be obtained by using Lemma [6] Here, use is also made of the antisymmetry of the 
CCR matrix 0. Now, application of Lemma [7] leads to sur(AjX) =( (AjX) 2 = AjXX T Ao, which, in 
combination with (l40l) . implies that 

ZZ T =<! -40A o A o t XX t A o A o t 0 = -40 AoA o t XX t A o A o t 0. 

Therefore, the operator inequality (l20l) is satisfied for d — 1 with F] = -^=0AoA ( } and Mo = 0. According 
to Theorem |TJ the perturbation Hamiltonian H\ being considered belongs to the uncertainty class Z/ 
with respect to which the system is robustly mean square stable, provided the following LMI 

A T n+ TL4 + (/ii + y)n H-110Ao 11 Ft Ao AfJ" =( 0 

Mi 

(obtained from (l22l) f holds for some n >- 0, Ml > 0, 7 > 0. 

Example 2: Consider a perturbation Hamiltonian H\ Y, k =\ r k cos(X^X + (j> k ), where X k G W are 
given vectors of spatial frequencies, r k > 0 are amplitudes and 0 ^ (j) k < 2n are initial phases which form 
complex amplitudes a k r k e l< ^ k . The corresponding Fourier transform in (J9]) is /i(A) = (<M<5(A — 

X k ) -\-a~k5 {A +A^)). Then the perturbation vector Z admits the decomposition (1341) with unit coefficients 
c k — 1, where, in view of Example [H 

Z, = ir k ®X k {^ X+M 

= -2r k ®h sin( A^X + (j) k ), k=l,...,d. 
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Hence, by applying Lemma [7] twice, it follows that 

Z k Zj = -4^04* sin 2 (A** + <fe)\ T © 
r2©4(4 T Z + ^) 2 4 T 0 

“4r^0A)t ^(1 + Qk'j^'k®’ 

where (O k are arbitrary positive real parameters. Therefore, the inequalities (031) are satisfied with the 
following parameters: 

^k'—^kxl * ®W, Ak) k'.=4rl(j) k --1 © 1 2 . (41) 

y Mi 

These can be employed in order to find the parameters T k and po according to (1361) and 071) of Lemma [8] 
and then proceed to the robust mean square stability analysis through the LMI (l22l) as described in 
Theorem [T] 

Example 3: Let E := H\ — L k =\ r k cos,(XjX + <j> k ) be an error of approximation of the perturbation 
Hamiltonian by a trigonometric polynomial from Example [2] Suppose its contribution i[E,X] to the 
perturbation vector Z satisfies 

-[E,x][E,x] T ^mrxx T r T +ni a , 

with r G R' 7X ” and p G R. A combination of Lemma [8] with the results of Example [2] leads to an 
augmented set of parameters which consists of <f>i,... ,d> f /, poi ,..., pod from (HTT) and <t>d+i := T, 
Po.d+i := P- The remaining part of the robust stability analysis procedure is carried out as before. 

X. CONCLUSION 

In this paper, we have presented a novel model for perturbations of Hamiltonians in a Weyl 
quantization form for a class of open quantum stochastic systems with linear coupling to the external 
boson fields. The time evolution of weighted mean square functionals of the system variables and a 
related dissipation inequality have been studied in order to develop sufficient conditions for robust mean 
square stability. An admissible class of Hamiltonian perturbations of a given stable linear quantum system 
has been formulated to guarantee stability of the resulting perturbed system. We have also discussed 
feasibility of these conditions in terms of the Weyl quantization model. This approach to the modelling 
and robust stability analysis of uncertain quantum stochastic systems has been demonstrated for several 
examples with Hamiltonian perturbations in the form of trigonometric polynomials of system variables. 
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